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Baker, Ball and Zachariasen have developed an elegant formulation of the dual superconducting 
model of quantum chromodynamics (QCD), which allows one to use the field equations to eliminate 
the gluon and Higgs degrees of freedom and thus to express the interaction between quarks as an 
effective potential. Carrying out an expansion in inverse powers of the constituent quark masses, 
these authors succeeded in identifying the central part, the spin-dependent part and the leading 
relativistic corrections to the central potential. The potential offers a good account of the energies 
and splittings of charmonium and the upsilon system. Since all of the flavor dependence of the 
interaction is presumed to enter through the constituent masses, it is possible to test the potential 
in other systems. Logical candidates are the heavy B-flavor charmed system and the heavy-light 
systems, which should be more sensitive to the relativistic corrections. Lattice gauge calculations 
furnish an additional point of contact for the components of the BBZ potential. Some preliminary 
calculations of the energies of B and D mesons are presented and the challenge of agreement with 
experiment is discussed. The spinless Salpeter equation is used to account for the effects of relativistic 
kinematics. 
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INTRODUCTION 

After an extensive investigation Baker, Ball and Zachariasen (BBZ) have developed a number of arguments 
to support the use of the dual color fields (derived from the color electric vector potential C^), instead of the usual 
Yang-Mills fields (derived from the magnetic vector potential A^), to formulate an effective theory of long-distance 
quantum chromodynamics (QCD). The Dyson-Schwinger equation may be used to justify some of their arguments, 
since it shows that corrections to the leading terms in the propagator A.4 are more singular in the infrared region 
that the leading term. This behavior is not true for the propagator, and thus one should use the dual fields in the 
long-distance region. Furthermore, BBZ argue that the connections between the dual fields and the Yang-Mills fields 
are nonlocal and require the introduction of an additional field. They use this requirement to good advantage by 
introducing Higgs fields, which allow a symmetry breaking mechanism to tmerate and thus confinement is explained 
as a Meissner effect, in the spirit of the ’t Hooft-Mandelstam conjecture |^,Q, where a magnetic condensate exerts 
pressure on the gluon field lines. The BBZ Lagrangian leads to an effective field theory with many of the properties 
expected from QCD, such as unitarity, renormalizability and non-Abelian gauge invariance. One can also show that 
the theory undergoes a deconfinement transition |^. A recent examination of the foundations of the BBZ model 
led to a reformulation of the theoretical underpinnings in a Wilson loop context. 

In their more recent work BBZ have constructed an interaction potential between a quark source and its 

antiquark by developing an expansion in inverse powers of the constituent masses. They use the field equations to 
eliminate the gluon and the Higgs fields and express the interaction energy in terms of the quark and antiquark 
variables. They can readily identify the central potential, the spin-dependent potentials and the leading relativistic 
corrections to the central potential. In this paper we explore some of the problems that arise when one attempts 
to use this potential in systems where relativistic corrections are more important. In particular, we attempt to give 
a simultaneous account of the properties of charmonium, the upsilon system and the heavy light systems using the 
spinless Salpeter equation for relativistic kinematics and a single set of potential parameters 00- In addressing 
this issue we will be led to examine questions about the relationships between the constants in the BBZ potential 
components and the parameters associated with mass renormalization [T^ . We note in passing the interesting work 
of Crater and Alstine @0 that explores the question of a universal quark-antiquark potential in the context of the 
Todorov equation and constraint dynamics. 

Since spin dependence is not a natural companion to work based on the spinless Salpeter equation, we will focus 
on spin-averaged energies as the point of contact with experiment. This goal is consistent with the goal of the 
BBZ investigations to produce an effective theory of long-range QCD, since spin-dependent splittings are presumably 
more sensitive to the short-range behavior of the potentials, where running coupling constant effects are important. 
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Although it is possible to accomodate such effects within the framework of the BBZ model, this involves additional 
assumptions. We have plans to pursue the role of the spin dependence in the context of the work of Hardekopf 
and Sucher |17[ . 

Our method of solution of the spinless Salpeter equation is based on a recent improvement |l^J^ of the Rayleigh- 
Ritz-Galerkin method. The purpose of our calculation is to see how successful one can be with the leading relativistic 
corrections and what role experiments might play in constraining these corrections and their relativistic generalizations. 
We establish contact with recent precision lattice gauge calculations of Bali, Schilling and Wachter |20 ^ in order to 
have an additional guide for evaluating the relativistic corrections. 

The highpoints of the BBZ formalism are summarized in Sec. ||, and the solution of the spinless Salpeter equation 
is presented in Sec. Some of the problems associated with the Salpeter wave functions and the salient features of 
the interpretation of the momentum operators are discussed in Sec. III. Our results and conclusions are presented in 
Sec. |f^ 


I. BBZ MODEL 


The fundamental degrees of freedom of the BBZ model are an octet of dual gluon potentials (7“ (o = 1,..., 8 ) and 

a triple of scalar Higgs fields Because of the symmetry breaking mechanism, these Higgs fields have nonzero 
vacuum expectation values, which give rise to effective masses for the dual gluons and two of the scalar Higgs fields. 
These effective masses are responsible for the exponential damping of all fields outside of the region of the sources 
and the flux tube that connects them. The model is based on the Lagrangian density, 


C = TR 


fH^ -D^ 


) + 


-W[B), 


( 1 ) 


where all the field variables in Eq. (|^ are scalar products with the SU(3) color matrices Aa (e.g. ^ E C^Xa) and 

W{B) denotes the Higgs potential. The covariant derivative in Eq. (j^). 


D^B = dfj,B - ig 


Cf.,B 


( 2 ) 


is necessary for non-Abelian gauge invariance. The electric displacement D and the magnetic field H are related to 
the vector potential by the non-Abelian relations, 


D = -V X C - 


ig 


[C,xC], 


(3a) 


H = -VCo-9oC-f5[C,Co]. (3b) 

To derive the simplest set of field equations with non-Abelian couplings, BBZ choose a gauge where the components 
of the Higgs triple point in three independent color directions, namely, B = (AyHi, —A 5 H 2 , A 2 B 3 ) and all the dual 
gluon field variables are proportional to the color matrix Y = X^/\/Z. This particular choice leads j2^ to finite results 
for the commutators of the first two components of B in Eq. (|^), but the remaining commutators in the equations 
above vanish. Thus all the relationships in the gluon sector are Abelian. Evaluating the trace of Eq. (0) yields 

£ = ^ (tf - D^) -H 2 d^B ■ d^^B + 2 {bI + B^) + W{B). (4) 

Since this Lagrangian density is symmetric in the components Bi and B 2 , we may choose them to be equal. The 
fourth-order term in Eq. (^ leads to effective mass terms in the field equations of the dual vector potential (7^ and 
the first two components of B, thus restricting these fields to the flux tube region and its vicinity. 

To complete the calculation of the interaction energy of a quark and its antiquark, one must introduce these particles 
as classical sources with masses mi and m 2 and spins (Ti and tT 2 , which are located at xi and X 2 and connected by 
a Dirac string. It is important to isolate the string contributions to the electric displacement and the magnetic field, 
Ds and Hs, which are produced by the sources, 

p(x) = e [(5^(x - xi) - (5^(x - X 2 )] T, (5a) 
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(5b) 


j(x) = e [vi<5^(x - xi) - V2 i5^(x - X 2 )] Y, 

where the velocities Vi are simply the time derivatives of the quark coordinates Xj. The relationships between the 
fields and the potentials also include explicit contributions from the color magnetization M and the color polarization 
M, which arise from the spin dependence of the sources. Thus, 

D = -V X C + Ds - P, (6a) 


H = -VCo - doC + Hs + M. (6b) 

From the Lagrangian density of Eq. (^) one can derive the field equations for the vector potential and the Higgs 
fields B. Their couplings to the sources are given by Eqs. (^, and thus the energy content of the fields surrounding 
the sources is completely determined by the solution to these field equations. Substituting these solutions into the 
expression for the Lagrangian density and integrating over all space gives the energy of the quark-antiquark system 
as a function of the source parameters r = Xi — X2, v; and cTj. Thus the potential is given by 


H(r,Vi,(Ti) 


J £(x, Vi, (Ti, C, i9jC, B). 


(7) 


It is straightforward to identify the central potential, the spin-dependent potential and the velocity-dependent 
corrections to the central potential, that is. 


V = Vo + VsD + 


( 8 ) 


where the central potential is given by 

Vo = Ar- - 0.646(9) 

3r 

It is not necessary to list the spin-dependent potential here, since it will not contribute to the spin-averaged energies 
in the lowest order. The velocity-dependent corrections to the central potential are given by the sum of an angular- 
momentum term and a radial-momentum term, namely, Vy 2 = Vang + Vrad- Each of these terms consists of two parts. 


Vang = V+ [r X (vi - V 2 )]^ /4r^ -h V_ [r X (vi -k V2)]^ /4r^, 


(10a) 


Vrad = V|| [r • (vi -k V2)]^ /4r^ -k Vl [r • (vi - V2)]^ /4r^. 


(10b) 


The potentials listed in Eqs. (10) are not all independent. Because of the Lorentz covariance of the theory ]9|,|13|, two 
of them may be expressed in terms of the central potential as follows: 


H- = 

BBZ 0 list explicit results for the other two, 


'll--if' 


14 = _^g-i.iVA/as r _ o.208Hr -k 
3r 


( 11 ) 


(12a) 


Vl = -k 0.0885(12b) 

dr 

Thus only two parameters, the string constant A and the strong coupling constant as are required to determine the 
potential. Of course, in the center-of-momentum frame the constituent masses become part of the parametrization 
through the relations, Vi = p/mi and V 2 = p/m 2 . Since agreement with experiment will not be easy to achieve, it 
is important to explore the sensitivity of the calculation to various parts of Vy 2 . Thus we report results below where 
the form of the longitudinal component is changed to 


vL = 


4as 

dr 


(13) 
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and to 


Vi 


" = _ 1^^-0.6851/17^ 


3r 


(14) 


In the limit of perturbative QCD {A 0) the potentials V+, V-, Vji, Vl all reduce to multiples of the Coulomb 
potential. Thus, the central potential and its leading relativistic correction reduce to the Darwin limit 


Vo + K2^-^ 
3r 


1 _ _ (vi . V2 + Vl • f V2 • f ) 


(15) 


We determine our parameters in the heavy quark systems, where the expansion in powers of the inverse masses is 
expected to have the most validity. Using the spin-averaged IS, 2S and IP energies of the upsilon system |24| and 
the IS energy of charmonium, we find that 


A = 0.235 GeV^; -as = 0.398; mt, = 4.710 GeV; = 1.320 GeV, 


(16) 


for the case that Vl is given by Eq. (H) above. 


II. SOLUTION OF THE SPINLESS SALPETER EQUATION 


Our semirelativistic potential model is based on the Salpeter equation, where 


H = 


mt 


-k 





(17) 


and the potential is given to O(v^) by Eq. (||). This model may be thought of as a three-dimensional reduction of the 
Bethe-Salpeter equation We will find the exact solution to this equation with the central potential of Eq. (|) 

using the approach of Refs. [ |18| , p^ to the Rayleigh-Ritz-Galerkin method, where the wave function 'I'n is expanded 
in terms of a complete set of basis functions. Then we will add the effects of Vy 2 with a first-order perturbation 
improvement. 

The basis functions |p^ , |l8| , |l^p^ used for the calculation include a spherical harmonic factor Y^ra (?) and the radial 
functions. 


Rni = iV„,/33/2(2/3r)^e-^"L^^+"(2/3r), (18) 

where = 8(n!)/r(n + 2^ + 2), /3 is a scale parameter, and the oscillations necessary for completeness are supplied 
by the associated Laguerre polynomials L^^+^(2/3r). Two important advantages accompany the use of these basis 
functions. The first of these is that all of the basic matrix elements that we require to find the eigenvalues and 
the eigenfunctions of Eq. (|l7|) can be generated from straightforward analytic expressions. The second is that the 
confinement factors of Eq.^8|) are a sufficienly good approximation to the actual wave function that the sizes of 
the requisite matrices are relatively modest (20 x 20 or 40 x 40). Analytic expressions for the matrix elements of 
the operators (p^)nn', the linear potential, (Ar)„„/, and the Coulomb potential (—4as/3r)„„/ are listed in Ref. [^ . 
There it is also explained how one can calculate the matrix elements of the nonlocal kinetic energy operators in Eq. 
(0) by effectively taking the square roots of the matrix representations of the operators 4- . 

Thus the remaining challenge in finding the solutions to the Hamiltonian operator of Eq. (f^) resides in the Yukawa 
factor present in the central potential of Eq. (^. We write the central potential as a Cornell Potential (linear + 
Coulomb) and a residual part by adding and subtracting a Coulomb piece, that is, 

Uo = Vcoruell + ^ (l " > (19) 

where ry = A/as- To determine the matrix elements of the residual part, we use the following algorithm: 

1. Generate the matrix elements of the Top, that is, the (r)nn'- 

2. Diagonalize this matrix to find the eigenvalues r^, which give a unique signature of Top in the chosen Hilbert 
space. 
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3. From the diagonal elements form th e q uantities f{rk) = (1 — e /r^, which constitute a signature of the 
residual part of the potential of Eq. 

4. Restore to the original basis with the unitary transformation that connects the original matrix representative 
of Top to its diagonal signature. 

5. Examine the procedure for stability by varying the sizes of the matrices. 

In order to have a definitive test for the stability of this five-step algorithm, it is convenient to have analytic results 
for some of the low-lying S-state matrix elements. With the basis functions of Eq. (|l^) it is straightforward to evaluate 
a few matrix elements of the residual operator f{r) when £ = 0. We list five of these from the first two rows of the 
matrix, 


,,, l36{2 + 5) ,,, /3<52(3 + <5) 

(1 + <5)2 ’ (/)oi V3(l + 5)3’ 

(20a) 

,,, £35^3 + 5) 

76(1+ J)4’ 

(20b) 

pSi2 + 3(5 + 4,52 + (53) 

0 ) 11 - (l-f(5)4 

(20c) 

,,, /352(4 + 4(5 + 552 + 53) 

72(1+ <5)3 

(20d) 


where S = r]/{2f3). We have verified that the matrix elements of Eq. (20) have the correct small ry behavior, where a 


power series expansion is valid, and the correct large ry behavior, where the matrix elements must reduce to those of 
the Coulomb potential. 

The results for the stability study of the matrix elements of the function /(r) are presented in Table |[ There results 
are listed for five different matrix sizes used to find the eigenvalues of Top, as discussed in step 2 of the algorithm 
below Eq. (HI)- Since there is no discernible difference in the matrix elements listed in the middle four columns and 
these results agree with the analytic results, one concludes that using a 40 x 40 matrix should generate very accurate 
results for {f)nn' ■ The energy eigenvalues listed in Table | were obtained with kinetic energy matrices the same size as 
the potential energy matrices. For the three smaller cases the Hamiltonian matrices were also chosen the same size as 
the kinetic and potential energy matrices, but for the two larger cases a 20 x 20 matrix was used for the Hamiltonian 
matrix. All of these results are consistent with our earlier study ||^ of the stability of the algorithm used to take the 
square roots of the kinetic energy operator and the choice of 20 x 20 as adequate for the Hamiltonian matrix. 

In summary our calculation of the energy eigenvalues and eigenfunctions requires four separate matrix diagonaliza- 
tions, two to determine the kinetic energy operators, one to determine the residual operator f{r), which represents 
the difference between the Yukawa potential and the Coulomb potential, and a fourth to diagonalize the Hamiltonian 
matrix. For the first three diagonalizations we use 40 x 40 matrices, and for the last we use a 20 x 20 matrix. One 
can generate accurate results for a range of scale values. The choice [3 = 2.0 GeV is a good one for the upsilon system, 
but (3 = 1.0 GeV is better for the heavy-light systems. 


III. VELOCITY-DEPENDENT CORRECTIONS 


Working in the center-of-momentum frame, where vi = p/mi and V 2 = —p/m 2 , allows us to write the velocity- 
dependent potentials of Eq. (10) in the form, 


Tanp — 


_L + _L)'IrA_L__L)'F_ 

rui Tri2 J 4r^ K^rii 'm 2 } 4r^ 


(r X p) • (r X p), 


(21a) 


^rad — 


1 

mi 


1 

m 2 


4j,2 


1 

mi 


1 

m 2 


4p2 


(r-p)(r-p). 


(21b) 
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It is important to note that Eq. (21) is a classical expression for the velocity-dependent potential energy of two 
classical sources, and thus it carries no instructions for ordering the momentum and position operators. Certainly, 
the quantum mechanical generalization of Eq. (H must respect the principle of Hermiticity. Although this principle 
is an adequate guide to determine the appropriate combination of any two noncommuting operators, some additional 
considerations are required for more complicated combinations. To resolve the a mbig uity, we choose to interprete the 
combinations of momentum and position operators in Eq. (|^) with the Gromes ||27| double-bracket notation. 


{{piF{r)pj}} = - [piFpj +pjFp^+piPjF + FpiPj]. 


( 22 ) 


This expression has the virtue that it is manifestly Hermitian and that it is involved in the nonrelativistic reductions 
of several covariant Dirac potentials. 

One additional difficulty requires some comment. The S-state expectation value of V), 2 , the first-order perturbation 
theory improvement, leads to singular integrals. Such terms arise because of the dominant influence of the Coulomb 
part of Vo upon the Salpeter wave function at short distances. Nickisch, Durand and Durand |2^ have shown that 
the small r behavior of the IS radial function is given by Ris{r) Br~P, where p « 0.25 for realistic values of as- 
To illustrate this behavior we choose one of the factors that arises from the last term of Vrad, that is 


II 



(23) 


The small r behavior of the integrand is given by G{r) which leads to a divergent integral at 

the lower limit. An important part of the reason for this divergence can be traced to the nonrelativistic limit used 
for the velocity operators. The relativistic expressions for the velocity operators of Eq. (|^ are Vi = pj\J‘ml + 
which have a finite large-momentum limit, in contrast to Vi « p/mi, which we used to obtain Eq. Thus, the 

constituent masses should provide a high-momentum cutoff in a more accurate treatment. Translating this observation 
into position space suggests that the Compton wavelengths of the constituent quarks should provide a short-distance 
cutoff, which would prevent the singular behavior. Since the singular behavior arises in the perturbative limit of 
Eq. and this equation is symmetric in the two masses, we choose a cutoff that is the geometric mean of the 
constituent Compton wavelengths, that is, A = (mim 2 )^^^. This choice does not introduce any new paramters into 
the calculation. Since all of the singular integrals that arise in the calculation of the expectation value of 1^2 have 
exactly the same behavior, one cutoff suffices to eliminate this behavior. 


IV. RESULTS AND DISCUSSION 


Our results for the IS B meson masses are shown in Fig. |^, where they are compared with the spin averages 
of the measured B meson masses [^. The full Salpeter curve shows the results of a calculation with the central 
potential of Eq. (||) only. This curve exhibits a uniform decrease, as the light quark masses decreases, in contrast 


to the pathological behavior |25 ^ of the Schrodinger equation results. It suggests that the Salpeter equation is a 


reasonable starting point for a potential model calculation since it is not too far above the data. The remaining four 
curves of Fig. include the effects of the velocity-dependent potential 1^2 of Eq. Each of these curves include 
perturbative and nonperturbative effects since they use the same expressions for V+, V- and Vj|. The dotted curve is 
based on the last of Eqs. (for Vl- Thus it includes all of the physics described by the BBZ model. In contrast, the 
dashed curve includes only perturbative effects in Vl, since it is based on Eq. The relativistic correction of the 

dotted curve has the wrong sign to move the Salpeter results towards agreement with experiment. The appearance of 
a pronounced minimum near 0.325 GeV shows the return of nonphysical behavior similar to the Schrodinger results, 
that is, the mass of the ^stem increases as the constituent mass decreases. Such behavior has also been noted by BBZ 
in an Appendix of Ref. |^. However, the results of the dashed curve exhibit the characteristics necessary for agreement 
with experiment. The dashed curve intersects the measured values near mu = 0.180 GeV and Wg = 0.220 GeV, which 
are reasonable values for the constituent masses . 

The remaining two curves in Fig. are based on additional modifications of Vl or 142 , which might be helpful 
in isolating the reasons behind the stri king differences between the dotted and dashed curves there. The dot-short 
dashed curve of Fig. || is based on Eq. ([l^. All of the difference between this curve and the dashed curve is due to 
the exponential damping of the Yukawa factor. Such a pronounced dependence is a consequence of the weighting 
of the integrand of the cutoff form of the integral of Eq. ( |^ ) . 

The long and short dashed curve is the result of an attempt to see if the freedom implicit in mass renormalization 
can be helpful in removing the discrepancy between the BBZ model (dotted curve) and experiment. Barchielli, 
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Brambilla and Prosper! [|^ show that a redefinition of the constituent masses leads to changes in the components 
of the velocity-dependent potential Vy 2 —at least in the nonrelativistic limit—, since, if one makes the replacement 
rui rrii + A, the Hamiltonian 

2 2 

H = mi+7712 + -1- - \-Vo + VsD + Vy 2 , (24) 

2777-1 2777-2 

is changed to 

HR=mi+m2 + 2A + i^ + ^-^ (^ + ^) + Vo + VsD + Vy2. (25) 

2toi 2 to 2 2 \mf m^J 

The two A-dependent terms can be absorbed in a redefinition of the potentials Vq and Vy 2 . However, one must 
be careful to respect the momentum dependence prescribed by the form of Eq. (^^. This requires that the same 
constant be subtracted from each of the component potentials, V+, V_, Vn and Vl- The potential parameters obtained 
by requiring the “mass-renormalized ” component potentials to fit the upsilon energy differences are identical to those 
obtained with the full BBZ potentials, as expected. The differences between the long and short dashed curve and 
the dotted curve, which become larger at smaller values of the light quark mass, probably reflect the fact that one 
requires a different mass renormalization procedure in this region. 

It is worth noting that the parameters used to obtain each of the curves in Fig. are slightly different, since in 
each case they were obtained by fitting the energy levels of the heavy quark systems. For example the dotted curve 
is based on the set, A = 0.229 GeV^, as = 0.307, mi, = 4.717 GeVand me = 1.308 GeV, which differ slightly from the 
values listed in Eq. (p^, the set used to calculate the dashed curve. The potential parameters for the two remaining 
curves in Fig. p] are identical to those for the dotted curve. 

The results for the IS D meson masses depicted in Fig. || are very similar to those of Fig. |^, since each curve there 
has the same qualitative features as its counterpart in Fig. |^. The dashed curve in Fig. || intersects the measured 
masses near m„ = 0.130 GeF and mg = 0.200 GeF. These values differ somewhat form those determined in Fig. ||, 
but such differences could easily be accounted for by a running mass. 

Results for the P state masses of the D mesons are presented in Fig. ||. The dotted curve has a shallow minimum 
around 0.275 GeV, in qualitative accord with the S state behavior. The dashed curve again intersects the data points, 
but the values that one obtains for the constituent masses, namely m„ = 0.110 GeF and TOs = 0.125 GeF are 
uncomfortably close together. The remaining two dashed curves in Fig. have qualitative differences from their S 
state counterparts. 

Thus, our conclusion is that we have not been able to use the full BBZ potentials to simultaneously account for the 
energies of the heavy quark systems and the heavy flavour systems with a single set of potential parameters in the 
context of the spinless Salpeter equation. It will be very interesting to see if calculations based upon the formalism 
of Hardekopf and Sucher 0 will remove this discrepancy. 

Since lattice gauge calculations are becoming more and more accurate, it is worthwhile to examine the recent 
calculations ^ of Bali, Schilling and Wachter (BSW) to see if they contain any guidance about which form of Vl is 
best to use. BSW present results for the velocity-dependent spin-independent corrections F,, Fc, Vd, Ve, which are 
defined as in the recent work of Brambilla and Vairo . These component potentials are certain linear combinations 
of the potentials F+, F_, V\\, Vl- For example, Fc and Fg are defined by the relations, 

Fc = (F_ + Fl - F+ - F|) /2, (26a) 


Fc = (F+ + F_ - k)| - Vl) /4. 


(26b) 


In Fig. U we compare two of our determinations of Fc with representative lattice data taken from Fig. 15 of Ref. 
plf . The dotted curve there is based on the use of Eq. ( p^ for Vl and the paramter set listed above in this section. 
The dashed curve there is based on the use of Eq. (p^ ) for Vl and the parameter set of Eq. (^). The lattice data 
does seem to favor the dotted curve over the dashed curve, but the error bars in the lattice data are large enough to 
prevent any strong indication. It is worth noting how much more sensitive the energy level calculations are to the 
differences between these two cases than the potential of Fig. 
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FIG. 1. Masses for IS B mesons as a function of the light quark constituent mass. 


FIG. 2. Masses for IS D mesons as a function of the light quark constituent mass. 


FIG. 3. Masses for IP D mesons as a function of the light quark constituent mass. 


FIG. 4. Comparison of the lattice results for 14 with those based on two different forms of Vl used in our calculations. 


TABLE I. Matrix elements of the residual potential function /(r) and S-state energy eigenvalues for the upsilon system. The 
parameters used for this calculation are those listed in Eq. (p^ and /3 = 2.0 GeV. 


Matrix elements 

5x5 

10 X 10 

20 X 20 

30 X 30 

40 X 40*^ 

Analytic results 

(/)oo 

0.38648 

0.38648 

0.38648 

0.38648 

0.38648 

0.38648 

(/)oi 

0.03346 

0.03346 

0.03346 

0.03346 

0.03346 

0.03346 

(/)02 

0.00318 

0.00318 

0.00318 

0.00318 

0.00318 

0.00318 

(/)03 

0.00031 

0.00031 

0.00031 

0.00031 

0.00031 


(/)ll 

0.35304 

0.35304 

0.35304 

0.35304 

0.35304 

0.35304 

(/)l2 

0.04799 

0.04799 

0.04799 

0.04799 

0.04799 

0.04799 

(/)33 

0.30244 

0.30264 

0.30264 

0.30264 

0.30264 


(/)44 

0.26930 

0.28336 

0.28336 

0.28336 

0.28336 


E{IS) 

9.47430 

9.47368 

9.47344 

9.47340 

9.47341 


E{2S) 

10.00877 

10.00718 

10.00703 

10.00704 

10.00701 


E{3S) 

10.38178 

10.36697 

10.36683 

10.36683 

10.36682 



'^The Hamiltonian matrix used for the eigenvalues in the last two columns was 20 x 20. 
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